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A  methodology  for  analysis  of  periodic  con¬ 
tent  in  sequentially  sampled  data  is  described. 
Whereas  previous  chronobiological  analyses  have 
been  limited  to  determining  global  parameters, 
multiple  complex  demodulation  (MCD)  provides 
local  properties  of  amplitude  and  phase  for 
each  datum  for  each  selected  demodulating  fre¬ 
quency  .  Real-time  occurrences  of  maxima  and 
minima,  exact  period  calculation  may  be  deter¬ 
mined  tron  rvmodu I.ites .  The  entire  analysis 
is  computerized  and  graphic  plots  are  output 
to  enhance  visualization  of  the  utility  of 
this  information. 


INTRODUCTION 

Mathematical  analyses  of  biological  and 
physiological  data  suspected  of  Inherent  periodi¬ 
cities  have  been  subject  to  various  methods  of 
time  series  analyses,  depending  on  the  nature  of 
the  information  desired.  Investigators  concerned 
vlth  phenomena  attributed  to  circadian  variations 
have  depended  greatly  on  the  Coslnor  method  deve¬ 
loped  by  Halberg  et.  al . *  Those  Involved  In  the 
detailed  examination  of  EEC  traces  have  fre- 

2  3 

quently  made  use  of  spectral  analysis  1  In  deter¬ 
mining  the  relationships  among  records  of  simul¬ 
taneous  measurements  from  several  electrode  pairs. 
Included  in  this  spectral  analysis  are  the  compu¬ 
tations  for  correlograma,  spectrograms,  perlodo- 

4  .  .  u  5,6, 7,8,9 

grams  ,  and  coherence.  Harmonic  or 

Fourier  analysis  have  been  used  by  still  others 
for  extracting  periodicities  not  Immediately 
apparent. 

Heart  rate,  EEC  activity,  eye  movements, 
wrist  and  arm  movements,  body  movements,  hormonal 
and  ensyme  secretions,  levels  of  biochemical 
constituents  are  stochastic  processes.  Their 
values  at  any  future  movement  cannot  be  exactly 
calculated  from  previously  known  measurements. 
Instead,  they  are  sampled  over  a  sufficiently  long 
Interval  to  detect  not  only  trends  and  abnormali¬ 
ties,  but  also  the  under lying  rhythms,  their  mag¬ 
nitudes  and  instantaneous  changes.  Then,  state¬ 
ments  may  be  made  as  to  the  probebf 11 ties  of  such 
events  occur lag  In  the  future. 

Harmonic  analysis  or  Interchangably,  Fourier 
analysis  is  applicable  where  It  la  aeetaed  that 
the  set  of  data  consists  of  a  series  of  sinusoidal 
components  irtilch  may  be  collectively  summed,  l.e., 
that  the  series  oscillate  more  or  less  harmonic ly 


over  time  according  to  these  components. 

Spectrum  analysis  provides  the  possibility  of 
any  particular  frequency  to  be  in  the  data  set, 
rather  than  setting  forth  exactly  the  actual  pre¬ 
sence  of  that  frequency.  Both  the  above  processes 
result  in  general  parameters  pertinent  to  the 
entire  ensemble. 

For  data  which  possess  gradually  changing 
amplitudes  and  phases  (hence,  implicit  changes  of 
frequencies),  applying  a  technique  that  restricts 
the  findings  to  only  a  single  frequency  and  its 
concomitant  amplitude  and  phase  clearly  limits 
the  evaluation  of  information  which  can  be  used 
to  enhance  interpretation  of  the  data. 

In  our  laboratory  we  have  utilized  extensive¬ 
ly  the  procedure  of  multiple  complex  demodulation 
(MCD)  in  analyzing  sets  of  data  as  diverse  as  EKG, 
body  temperature,  activity  movements,  EOG,  verbal 
communication  in  simulated  battle  environment,  and 
hormonal  secretion. 

The  complex  demodulation  technique  has  been 
set  forth  with  great  clarity  by  Bingham  et.  al. 
and  Bloomfield*1;  has  been  briefly  illustrated  by 
Walter1^,  applied  by  Childers*^,  Orr  &  Hoffman*^, 

and  more  recently  by  DePrlns  et.  al.1^.  There  is 
adequate  theoretical  basis  for  the  use  of  this 
method  in  extracting  from  a  sequential  data  series 
the  local  properties  from  each  epoch  in  the  time 
domain.  The  implication  Is  that  a  non-statlonary 
time  series  cannot  accurately  be  described  as 
having  parameters  which  are  constant  over  the  en¬ 
tire  set  of  data.  The  price  to  be  paid  for  this 

flexibility  according  to  Bloomfield11  is  less 
precision  in  the  description  of  pure  frequencies. 

METHOD  OF  MCD  ANALYSIS 

In  describing  each  of  the  steps  Involved  in 
the  multiple  complex  demodulation  technique,  the 
circadian  period  (24  hours)  will  be  used  in  the 

example. 

A.  The  Data  Set 

Measurements  sampled  sequentially  In  equal 
time  Increments  preferebly  for  at  least  2  cycles 
(48  hours)  will  comprise  the  date  set  depleted  as: 

X(t)  -  x j ,  x2 . Xj|  (1) 

where  each  x  element  la  an  epoch  value. 
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The  time  increment  may  be  one  minute  as  in 
EKG  data;  10  minutes  as  in  actigraph  (wrist  and 

arm  activity)  records1^;  or  integral  hours  as  in 
body  temperatures  and  hormonal  secretions.  The 
choice  of  sampling  interval  depends  not  only  on 
the  limitations  of  the  measuring  instrument,  but 
also  on  the  number  of  constituent  frequencies 
desired  for  analysis. 

Discrete  equal  sampling  intervals  introduce 

a  problem  known  as  "aliasing" , ^  in  which  fre¬ 
quencies  higher  than  the  Nyquist  frequency  are 
enveloped  by  lower  frequencies  with  which  they 
coincide  at  crossover  points.  The  highest  fre¬ 
quency  analyzable,  in  a  record,  then,  is  the 
Nyquist  frequency.  For  a  circadian  cycle  this 
limit  is:  24/2At,  where  24  represents  1  cycle 
in  hours,  and  At  is  the  sampling  interval  in 
the  same  unit. 


D.  Conversion  tQ  the  Complex  Domain 

The  first  step  in  the  actual  process  of  de¬ 
composing  the  data  set  into  its  constituent  fre¬ 
quencies  involves  the  conversion  of  the  values  to 
the  complex  domain.  This  is  necessary  for  the 
computation  later  of  the  proper  amplitudes  and 
phases.  In  the  conversion,  each  epoch  value  is 
multiplied  with  separately  (1)  a  cosine  wave  and 
(2)  a  sine  wave  containing  the  frequency  to  be 
elicited  (the  demodulated  frequency)  in  their 
arguments.  The  result  is  the  formation  of  2  terms 
for  each  data  epoch  as  follows: 

Z (Re)  *  7  *  Cos  Xt  (5) 

ri 

Z(Im)  -  7  *  Sin  Xt.  (6) 

^  j  * 


B.  Detrending  of  Data  Set 

The  data  set  is  detrended  by  subtracting 
the  mean  value  of  the  entire  set  from  each 
epoch  datum.  This  mean  value  is  commonly  refer¬ 
red  to  as  the  "zero  frequency"  and  the  result  is 
a  new  data  set  which  oscillates  sinusoidal-like 
around  a  zero  level.  This  set  is  now: 

Y(t)  *  yr  y2>  . yN  (2) 

where 

«- 1 . N  <3> 

t*l  u 

C.  "Folding-Out"  of  Data  Set 

Necessarily,  in  the  actual  analysis,  we  must 
deal  with  a  finite  set  of  time  series  data.  How¬ 
ever,  abrupt  truncation  of  an  infinite  series  on 
either  end  has  consequences  of  "end  effects" 
which  result  in  distortion  of  their  local  pro¬ 
perties.  To  minimize  such  effects,  we  have  ex¬ 
tended  the  length  of  our  data  set  by  "folding-out" 
the  epoch  values  at  both  the  beginning  and  ter¬ 
mination  of  the  data  series.  This  is  reasonable 

(Bloomfield11)  in  light  of  other  alternatives 
which  add  zeroes  to  either  side.  This  modifica¬ 
tion  results  in  yet  a  new  set  of  data  values: 

z(f0)  ■  v  Vi . v  Xi»  y2 . xM. 

V  yN-l . y, m-J  (4) 

where  m  is  the  number  of  folded-out  data  points 
and  f  -  o-l.  Note  that  the  underlined, 
are  the  beginning  and  end  respectively  of  the  data 
set  prior  to  the  folding-out  process. 

The  amount  of  data  points  folded-out  varies 
according  to  the  length  of  the  data  set.  Our  rule 
has  been  to  use  20X  of  the  total  for  a  long  data 
set  (  >  100  values)  and  5Z  for  short  segments. 

The  folded-out  data  set  becomes  the  working 
set  for  all  subsequent  mathematical  manipulations. 
The  final  output  however  retains  only  the  para¬ 
meters  of  the  original  epochs  for  statistical 
analysis  and  display. 


where  X  “  2irfj/S.R.#  f  ■  frequency  selected  for 

demodulation,  J  -  1,2,  ...,  k,  k  <_  S.R/2,  and 
S.R.  is  defined  as  the  saminling  rate  or  number 
of  sample  epochs  in  the  chosen  period.  For 
example,  if  the  sampling  interval  is  15  minutes 
and  the  frequency  of  2  cycles  in  the  basic  period 
(24  hours)  is  demodulated  then  X  -  2tt2/96. 

Implicit  in  each  datum  is  the  collective 
value  of  all  inherent  frequencies  contained  in 
the  data  set.  The  multiplication  by  the  sine  and 
cosine  terms  preserves  not  only  the  original  fre¬ 
quencies,  but  also  generates  sum  and  difference 
frequencies  between  the  modulating  frequency  and 
all  the  original  frequencies. 

In  the  Fourier  spectrum,  this  has  the  effect 
also  of  transposing  the  desired  modulating  fre¬ 
quency  to  the  zero  frequency  position.  This  will 
be  crucial  as  we  will  observe  later  in  the  filter¬ 
ing  process. 

Although  we  have  fixed  the  period  at  24  hours, 
the  multiple  complex  demodulation  allows  the  fre¬ 
quency  within  that  period  to  vary  and  so  we  can 
demodulate  each  frequency  of  interest  in  a  step¬ 
wise  manner. 

E.  Filtering 

The  extraction  of  the  desired  frequency  at  the 
zero  frequency  position  necessitates  the  exclusion 
of  other  components  near  to  it  in  the  energy  spec¬ 
trum.  These  may  be  noise  or  the  difference  fre¬ 
quencies  described  above  which  result  in  frequen¬ 
cies  near  the  demodulating  frequency. 

For  this  purpose,  an  exponential  filter  is 
used  and  is  passed  sequentially  through  all  the 
data  points  first  in  a  forward  pass,  then  a  re¬ 
verse  pass  and  the  whole  process  repeated.  The 
forward  pass  causes  a  shift  of  one  In  the  data  set, 
which  is  corrected  by  the  reverse  pass  thereby  pre¬ 
serving  the  true  phase  values. 

The  low-pass  characteristics  of  this  filter 
allows  the  power  at  and  near  the  zero  frequency 
In  the  spectrum  to  pass  through  while  excluding 
the  rest  of  the  spectrum.  Inevitably,  there  will 
be  some  "leakage"  of  power  through  the  filter,  in 
particular  from  whatever  residual  power  of  fre¬ 
quencies  located  adjacent  to  the  zero  frequency. 
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For  this  reason.  In  analysis  of  data  for  circadian 
rhythmicity,  the  epoch  values  obtained  from  the 
circadian  demodulation  and  properly  transformed 
into  remodulates  (as  will  be  discussed  later)  are 
subtracted  from  their  corresponding  values  In  the 
folded-out  data  set.  The  demodulations  of  all 
frequencies  greater  than  one  are  then  performed 
on  this  difference  set. 

The  filter  employed  in  our  procedure  con¬ 
sists  of  two  parts: 

F,  *  [  (U)  2  +  (V)2]*1 

F2  -  exp (a) , 

where 

(U)  -  1.0  -  exp (a)  *  Cos  (GF  •  T) 

(V)  -  exp (a)  •  Sin  (GF  •  T) 


The  remodulate  values,  after  truncation  of  the 
folded-out  epochs  at  the  beginning  and  termination 
of  the  data  set,  comprise  a  smoothed  function  of  the 
desired  demodulated  frequency  with  the  proper 
phases  intact.  Peak  and  trough  amplitudes  and 
their  corresponding  real  time  values  may  then  be 
determined  over  the  entire  length  of  the  data  series 
for  each  frequency  demodulated. 

For  a  record  consisting  of  several  cycles,  any 
shifts  of  trough  and  peak  occurrences  In  real  time 
may  be  detected  on  a  daily  basis.  In  addition,  the 

(7)  length  (in  hours)  of  the  actual  circadian  period 
may  be  calculated  from  peak  to  peak,  trough  to 

(8)  trough,  or  from  zero  cross  over  points  from  the 
remodulate.  Changes  in  period  lengths,  thus  also 
signal  phase  changes  and  may  be  confirmed  from  the 
phase  plot.  Also  to  be  noted  is  the  asymmetry 
between  the  positive  portion  of  the  cycle  and  the 
negative  portion. 


and 

T  =  2tt/S.R. 
a  -  -T 

GF  *=  variable 

and  SR  is  defined  as  in  (6).  The  gain  factor,  GF, 
may  take  values  from  0.1  to  0.9  depending  on  the 
magnitude  of  the  original  data,  i.e.,  the  larger 
the  data  values,  the  smaller  the  gain  factor  re¬ 
quired.  It  goes  without  saying  that  when  differ¬ 
ent  records  are  analyzed  for  direct  comparisons, 
the  same  gain  factor  should  be  used  in  all  cases. 

In  actual  practice,  the  filter  operation  in 
a  single  pass  involves  summing  a  proportion  of 
each  epoch  value  with  a  proportion  of  the  pre¬ 
vious  adjacent  epoch  value  in  the  following  way: 

0^  (Re)  -  (Re)  •  Fx  -f  •  F2  (9) 

0  (Im)  -  Z  (Im)  *  F.  +  Z  ,  (Im)  •  F0  (10) 
n  n  i  n—  i  2 

where  n  -  m,  m-1,  ...1,  ...N,  N-l,  . . . ,N -£  and  m 
and  t  have  bean  defined  in  (4).  The  outputs  from 
each  filter  pass  are  used  as  new  Inputs  for  the 
next  pass  going  in  the  reverse  direction. 


G.  Selection  Criteria  for  “True"  Frequencies 
Inclusion 

If  biological  and  physiological  rhythms  are 
expressions  of  the  energy  state  of  the  system,  then 
consideration  must  be  given  as  to  which  frequencies 
are  "true11  contributors  and  which  are  merely  noise 
components.  We  have,  in  our  laboratory,  establish¬ 
ed  rather  strict  criteria  for  the  inclusion  of  any 
particular  frequency  as  being  a  constituent  of  the 
total  energy.  Several  steps  are  taken  in  this  de¬ 
termination  process  and  they  are  listed  in  the  fol¬ 
lowing  order: 

(1)  The  remodulate  values  of  every  fre¬ 
quency  demodulated  are  scanned  for  their  peak  oc¬ 
curence  values  which  are  stored  in  arrays.  The 
maximum  value  of  all  these  peak  values  is  selected 
as  the  basis  for  scaling  into  20  grey  levels,  i.e., 
scale  factor  =  20. /maximum  value. 

(2)  The  peaks  of  each  remodulated  fre¬ 
quency  is  multiplied  with  the  scale  factor  and  the 
products  binned  into  the  appropriate  grey  levels. 
Products  less  than  0.5  are  binned  into  a  zero  level. 

(3)  The  total  population  is  determined 
from  the  sum  of  all  the  counts  In  the  21  bins. 


F •  Power,  Amplitude.  Phase .  and  Remodulate 

The  final  outputs  from  tha  filtering  opera¬ 
tions  may  now  be  used  for  computing  the  parameters 
or  local  properties  of  each  epoch.  These  are: 

Power:  P_  -  2.0  [0  (Re)2  +  0  (Im)2]  (11) 

n  n  n 

Amp:  A  -  2.0  [0  (Re)2  +  0  (Im)2]*4  (12) 

n  n  n 

Phase:  <6 

ARCTAN  S©n  (lm)/0n  (Re)]  (13) 

Reoodultte:  R  -  2.0  [0  (Im)  •  Sin  (X.t)  + 

n  n  t 

0fl  (Re)  •  Cos  aft>]  (H) 

where  the  unit  time  Increment,  t  «  0,1,...,  N+2n, 
and  f  Is  the  demodulated  frequency. 


(4)  Beginning  with  the  zero  level,  a 
cumulative  tabulation  is  initiated  with  determina¬ 
tion  of  the  percent  of  total  population  at  the  end 
of  each  level’s  cumulation. 

(5)  When  cumulation  at  the  level  contain¬ 
ing  90Z  of  the  total  population  is  attained,  then 
that  level  is  chosen  as  the  cut-off  value  for  deter¬ 
mining  inclusion  of  a  frequency  as  a  real  component 
of  the  data  set, 

(6)  Those  frequencies  possessing  a  peak 
value  in  their  remodulates  at  or  exceeding  the  grey 
level  cut-off  value  are  selected  out  as  part  of  the 
biological  system  and  are  not  considered  random 
noise. 

(7)  The  energy  spectrum  for  biological 
systems  may  be  considered  partitioned  as  shown  in 
Table  I.  On  this  basis  the  frequency  bands  will 
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be  synthesized  as  the  sum  of  remodulates  of  fre¬ 
quencies  selected  by  the  grey  level  process. 

Peaks  and  troughs  and  their  respective  positions 
In  real  time  are  also  determined  and  recorded. 
Scanning  down  such  lists  allows  immediate  detec¬ 
tion  of  the  location  (time)  of  unusual  events. 

(8)  A  graphic  presentation  of  the  grey 
level  values  above  the  cut-off  point  in  a  Time  vs 
Frequency  plot  can  illustrate  with  startling 
clarity,  the  wave-shaping  properties  of  groups 
of  frequencies. 

STATISTICS 

Since  the  multiple  steps  taken  in  the  demodu¬ 
lation  process  resets  the  values  of  the  output 
series  with  no  direct  or  calculable  relationship 
to  the  original  data,  it  is  not  possible  to  make 
a  direct  comparison  of  these.  But  since  in  a 
time  series,  the  important  considerations  are 
the  changing  nature  of  the  values  over  the  whole 
sample  series  and  not  the  values  per  se,  tb - 
statistical  comparisons  may  be  made  on  "norma¬ 
lized”  data  sets  to  capture  any  inter-relation¬ 
ships  in  synchrony  over  time.  For  this  purpose, 
we  have  scaled  each  epoch  value  In  each  series 
with  a  scaling  factor  obtained  by  the  ratio  of 
20:  difference  between  maximum  and  minimum  value 
in  that  series.  Accordingly,  the  following 
statistical  parameters  are  computed: 

(1)  Variance  of  the  original  data  set. 

(2)  Variance  of  the  circadian  fre¬ 
quency  and  of  each  frequency  band  synthesized. 

(3)  Co-variance  of  the  original  data 
set  with  the  circadian  frequency  and  with  each  of 
the  synthesized  frequency  bands. 

2 

(4)  R  -  Correlation  of  original  data 
set  with  the  circadian  frequency  and  with  each 
synthesized  frequency  band. 

(5)  R  -  Correlation  coefficients  of  (A). 

(6)  RPER  -  a  percent  value,  related  to 
a  goodness  of  fit  measure.  It  is  the  ratio  of  the 
sum  of  the  squares  of  differences  to  the  sum  of 
squares  of  maximum  differences,  where  each  maximum 
difference  value  is  equal  to  20.  A  value  of  I00£ 
would  indicate  no  fit. 

REFINEMENTS  -  SPLINE  FIT 

Peak  and  trough  selections  In  the  above  pro¬ 
cedure  will  always  be  limited  to  the  time  of  a 
sampled  epoch  whether  or  not  the  peak  occurs  there 
in  reality  or  between  sampling  Intervals.  Since 
it  is  not  known  where  the  true  peak  occurs,  we  have 

18 

applied  a  spline  fit  technique  to  the  reaodulate 
values  in  an  attempt  to  provide  s  better  estimate 
of  peak  location.  The  result  of  the  spline  fit  on 
a  data  series  is  a  smoothed  function  which  passes 
through  all  date  points  used  in  the  fitting  rou¬ 
tine. 


For  example,  if  the  sampled  data  consist  of 
epochs  taken  at  15  minute  Intervals  and  it  is  de¬ 
sired  to  locate  the  actual  peaking  at  better  than 
the  15  minute  interval,  then  a  spline  fit  around 
the  output  peak  (obtained  in  the  procedure  given 
above)  is  applied.  The  equal  interval  time  axis 
is  converted  from  unit  incremental  steps  to  the 
real  time  scale  —  in  this  case  to  units  of 
minutes.  Two  data  points  on  each  side  of  the 
peak  and  the  peak  point  (five  values  in  all)  are 
subjected  to  the  spline  fit  procedure  to  obtain 
the  parameters  for  the  smoothed  function.  These 
parameters  are  then  used  to  determine  the  value 
of  the  function  between  time  equivalent  epochs  on 
either  side  of  the  peak.  Hence,  we  can  now  deter¬ 
mine  values  of  the  function  at  one  minute  incre¬ 
menting  intervals  Instead  of  15  minutes.  The 
highest  value  computed  from  this  function  around 
the  peak  is  taken  to  be  the  best  estimate  of  the 
peaking  and  its  corresponding  time  as  the  actual 
time  of  peaking. 

Alternatively,  for  even  coarser  sampling  in¬ 
tervals,  i.e.,  every  A  hours,  a  spline  fit  may  be 
made  on  the  remodulate  values  and  the  values  of 
the  new  function  at  fine  resolutions  of  15  minute 
interval  increments  may  be  calculated.  The  re¬ 
sisting  values  comprising  a  "stretched"  remodu- 
lat«  series  may  then  be  scanned  for  their  peaks 
and  respective  times  in  15  min  intervale  instead 
of  4  hours. 

The  advantage  of  applying  the  spline  fit 
technique  is  that  we  are  not  locked  into  describ¬ 
ing  peaks  and  troughs  only  at  the  sampling  epoch 
time,  but  can  obtain  a  finer  resolution  of  peak 
time.  The  disadvantage  is  the  slight  distortion 
in  amplitude  of  the  peaks  and  troughs  as  a  conse¬ 
quence  of  the  requirement  that  the  function  must 
pass  through  each  data  point.  This  is  a  small 
price  to  be  paid  for  the  ability  of  exactly  locat¬ 
ing  critical  values  in  an  oscillating  system. 

COMPUTER  PROGRAMS 

Computer  programs  are  available  for  the  meth¬ 
odology  described,  in  which  the  input  required 
are  the  sequential  data  set  and  the  variables  of 
sampling  rate,  epoch  length,  g*.in  factor,  starting 
time  of  the  series,  and  number  of  frequencies 
(beginning  with  1  cycle)  to  be  demodulated.  The 
output  provides  values  of  mean  power  for  each 
cycle  so  that  a  "pseudo"  power  spectrum  may  be 
plotted.  Additionally,  all  the  maxima  and  minima 
of  each  remodulate,  their  real-time  occurrences 
and  amplitudes,  and  the  period  lengths  are  listed. 

Plotting  routines  are  used  to  display  the  re¬ 
sults  of  the  MCD  technique.  They  illustrate  the 
synthesis  of  the  remodulstes  Into  characteristic 
bands  of  physiological  significance  and  their 
summations  are  compared  with  the  original  data. 

The  major  modulator  (circadian  rhythm  with  period 
approximately  2A  hr)  is  drawn  superimposed  over 
the  entire  data  series  and  the  corresponding 
phase  Is  separstely  traced  beneath  it.  Finally, 
rest  and  activity  cycles  are  plotted  on  a  dally 
basis  to  show  any  incipient  phase  shifting. 

In  the  illustrations  to  follow,  the  sampling 
rate  for  a  24  hour  period  was  used  and  the  fre¬ 
quencies  arc  cycles  within  this  period. 
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FREQUENCY  (CYCLES/24  HOURS) 


Table  I.  Band  Partition  of  Frequencies  and  Kqulvalent  Periods 


Bands  Frequency 

Cyclea/Day 


period  Bands 

Hours  Hr/mln/sec 


Frequency  Period 

Cy  ices/Day  Hours*  hr/niin/see 


Circadian 

l 

24.0 

24/00/00 

Supra-Ultradian  21 

1.14 

1/8/34 

22 

1.09 

1/5/27 

Intermediate 

2 

12. C 

12/00/00 

23 

1.04 

1/2/36 

3 

8.0 

8/00/00 

24 

1.00 

1/00/00 

4 

6.0 

6/00/00 

25 

0.96 

0/57/36 

5 

4.8 

4/48/00 

26 

0.92 

0/55/23 

27 

0.89 

0/53/20 

Sub-Ultradian 

6 

4.0 

4/00/00 

28 

0.86 

0/51/26 

7 

3.4 

3/25/43 

29 

0.83 

0/49/39 

0 

3.0 

3/00/00 

9 

2.7 

2/40/00 

Very  High  30 

0.80 

0/48/00 

10 

2.4 

2/24/00 

31 

0.77 

0/46/27 

11 

2.2 

2/10/54 

32 

0.75 

0/45/00 

33 

0.72 

0/43/38 

U1 trad lan 

12 

2.0 

2/00/00 

34 

0.70 

0/42/21 

13 

1.8 

1/50/46 

35 

0.68 

0/41/08 

U 

1.7 

1/42/51 

36 

0.66 

0/40/00 

15 

1.6 

1/36/00 

* 

16 

1.5 

1/30/00 

Rounded  to  tenths  or  hundredths 

17 

1.4 

1/24/42 

18 

1.3 

1/20/00 

19 

1.26 

1/15/47 

20 

1.20 

1/12/00 

-4 

i 


"Zo 


f 


// 


&  5^ 


>- 


GREY  LEVEL  CUTOFF  •  3 


a  33*  1  1***3  3*3 
3)3  3  4  4)  «l«  414  3443  33 
3  3333  3 


3  33  3)  3)  3)3333 


3)  3443)443  3)3 

33333 


3)3)33)33 
334444443  33)33 

3  334)3  33  4  43  33)3  3)33  3 


))  4)|  41  11*4  MMII  4t  44  4)  4  )3  4)  4  34  3)  3  3  3 


3  3  3  3  3  3  3 

)))IS33)3)3) 

3)343)  3333 
3  3  3 


333)444443)3 
333)3)4*44)  3)33)  33) 


3  3  3 

I  3  3)3 

3  3  3 


3  3  3  3  3 
3)411144 
3)14)3 


3  )))  3  )  3  )  33434)3) 

III))  )!)))) 

)))  3)441144 


U  *  *  3  4  4 

10  V  )  )  3  3 

•  4  3  ) 


3  3  1)3 

)  •  •  3  •  f  I 

I  3  I  I  I  3  I  I  I 

3  4  3  3  1  • 

3  * 


3  3  3 

3  )  3 


•\  * 
B  * 

«  P 

:r 


1919  0119  0719  1919  IMS  0119  0719  1919  1919  0119  0719  1919  1919  0119  0713  1319  1919  0119  0719 

9219  0«  19  1019  1919  2219  0419  1019  1919  2219  0419  1019  1919  2219  0419  1019  1919  2219  04  19 


TIME  (24  -  HOUR  NOTATION) 


Fig.  1.  Rsstodulate  Maxima  of  dwodulated  frequencies  distributed  aa  grey  levels.  The  highest 
amplitude  Is  Indicated  by  20.  Ths  lowest,  3,  Is  ths  selection  criterion  level  of  significant 
components  for  Inclusion  In  the  10Z  or  lees  of  the  maxima  population.  The  original  data  are 
15  sin  epoch  recordings  of  vrlet  end  ere  movements. 
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Fig.  4.  Circadian  remodulata  superimposed  over  original  data  with  corresponding  phase  plot  and  time 
of  maxima  and  minima  occurrences.  From  the  wrist  and  arm  movement  records  of  a  shift  worker. 
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HI  activity 
CZU  PEST 


-'■«  <  M  *  90  1  00  .  0  00  (*. 90  H  OC  II  ID  It  SO  30  00  i:  C3  H  OC 


TIME  OF  OAV  *iTh  24  00  HP  NOTATION 


Fig.  5.  Rest/activity  cycles  determined  from 
zero  cross-over  points  of  circadian  rhythm 
in  Fig.  4.  Note  phase  shift  as  the  work 
shift  is  changed  on  day  18. 
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